High frequency dielectric and magnetic anomaly at the phase transition in NaV 2 Os 
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We found anomalies in the temperature dependence of dielectric and magnetic susceptibility of 
NaV2 05 in the microwave and far infrared frequency ranges. The anomalies occur at the phase 
transition temperature T c , at which the spin gap opens. The real parts of the dielectric constants e a 
and e c decrease below T c . The decrease of e a (except for the narrow region close to T c ) is proportional 
to the intensity of the x-ray reflection appearing at T c . The dielectric constant anomaly can be 
explained by the zigzag charge ordering in the ab-plane appearing below T c . The anomaly of the 
microwave magnetic losses is probably related to the coupling between the spin and charge degrees 
of freedom in vanadium ladders. 



PACS numbers: 64.70.Kb 75.10.Jm 77.22Ch 75.40. Cx 



I. INTRODUCTION 

Two inorganic compounds, CuGeC>3 and NaV^Os, 
were extensively studied as quasi-one-dimensional (ID) 
magnets containing the chains of spins S ~ h coupled by 
the antiferromagnetic Heisenberg exchange. Both mate- 
rials show ID behavior of the susceptibility at high tem- 
peratures and the phase transition into the spin-gap state 
at low temperatures. The lattice deformation appears si- 
multaneously with the spin-gap opening at the transition 
temperature T c . The lattice deformation includes the 
doubling of the lattice period in the spin-chain direction 
and the spin gap is believed to be a result of the corre- 
sponding alternation of the exchange constants along the 
chains. 

The phase transition in CuGe03 is generally consid- 
ered to be a spin-Peierls transition. The driving force 
of this transition is the collective spin-lattice instability 
with the gain in the energy oLthe spin chains exceeding 
the loss in the lattice energyLra The experimental data 
on CuGe03 are in a good agreement with the theoretical 
results for the antiferromagnetic spin- ^ chains coupled to 
the lattice.13 Some deviations from the ideal spin-Peierls 
behavior found for CuGeOa are also well-understood. 
They are related to a considerable interchain interaction 
(the intcijchain exchange being about 0.1 of the intrachain 
exchanges) and a strong next-nearest-neighhor exchange 
(about 0.3 of the nearest-neighbor exchange 



For some time NaV2 0s was considered as the second 
inorganic spin-Peierls system. The ID-magnetic struc- 
ture of NaV2 05 was first associated with the chains of 
V 4+ -ions (spin S—l/2) along the b-axis of the orthorom- 
bic cxistal separated by the chains of nonmagnetic V 5+ 
ionsBa This structure with inequivalent vanadium sites 
was recently questioned in Refs. |^-|lC|. On the basis of the 
new structure refinement studies of the high-temperature 
phase of NaV2 05 the new ladder.- type structure with 
equivalent V-sites was suggested. LrB The ladders are ori- 
ented along the &-axes with the rungs along the a-axis. 
There is one electron per rung and the spin-i chains are 
formed by the electrons localized on V-O-V molecular 
orbitals on rungs and coupled by the exchange interac- 
tion along the 6-dircction. In this structure the charge of 
V-ions fluctuates, its average .izalue being 4.5. 

Recent experimental data,EJ as well as the theoreti- 
cal arguments, EL3~£il suggest that the phase transition in 
NaV205 can be a result of some charge ordering rather 
than the spin-Peierls instability. Two kinds of such an or- 
dering were considered. In the model of Ref. the elec- 
trons below T c become localized on one leg of each laddex, 
as in the initially proposed high temperature structured. 
In this model one has to invoke an extra mechanism to ac- 
count iQE_the spin gap at low temperatures. In the second 
modeEj'E.3 the electrons (or the V 4+ -ions) form zigzags 
on each ladder. The opening of the spin gap in this sce- 
nario is a direct consequence of the charge ordering. The 
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optical spectra of NaV205E-3o were first interpreted as 
an indication on the presence of inequivalent V-sites even 
at room temperature. However, these spectra could also 
be explained by the presence of only a short range charge 
order with relatively slow charge fluctuations. An extra 
confirmation that the phase transition in NaV2C>5 is not 
an ordinary spin-Peierls transition comes from the study 
of thermal conductivity: in contrast to the spin-Peierls 
system CuGe03, the thermal conductivity of NaV2C>5 
has a huge anomaly below T c , which can be naturally 
explained by the charge ordering.EZI 

Further information about the charge subsystem close 
to the phase transition point can be gained by studying 
the dielectric constant. Measurements of the dielectric 
constant at the frequency of 1 KHz along three princi- 
pal directions were performed in Ref. An anomaly 
with a strong anisotropy with respect to the electric field 
orientation signifying the rearrangement of the charge in 
NaV2C>5 has been found. 

In this paper we present combined measurements of 
the real and imaginary part of the dielectric and magnetic 
susceptibility of NaX^Os at the microwave frequency of 
36 GHz and of the refractive index in the far infrared 
range. For comparison the measurements of the same pa- 
rameters were performed also for CuGeC>3. Our results 
on NaV2 05 confirm the presence of strong anomalies in 
the dielectric constant along the a- and c-directions, al- 
though the type of the anomaly for the c-direction is dif- 
ferent from the one found in Ref. [l^. The anomaly in e a 
is naturally explained by the picture of the zigzag charge 
ordering suggested in Ref. |l3], 14| Possible reasons for the 
anomaly in e c are discussed and further experiments to 
check this picture are suggested. 



II. EXPERIMENT 

We studied the complex microwave dielectric and mag- 
netic susceptibilities by measuring the temperature de- 
pendence of the resonance frequency and of the Q-factor 
of the microwave TE104 cavity of the size 20x7. 2x3. 4mm. 
The resonance frequency / of the empty cavity is 36 GHz. 
The presence of a sample shifts the resonance frequency 
of the cavity. When the sample in the form of a thin 
plate is placed at the maximum of the microwave electric 
field ~E m w, with the plane of the plate parallel to the elec- 
tric field, the resonance frequency of the cavity should be 
shifted by the value Sf: 



J ~ v ■ 



(1) 



Here v, V are the volumes of the sample and of the cavity 
respectively, e' is the real part of the dielectric constant. 

When the sample plane is oriented perpendicularly to 
the field E mu) , the shift of the frequency is given by the 
relation: 
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In this case the frequency shift is smaller due to the 
depolarization factor effect. Formulae (1) and (2) are 
valid in the quasistatic approximation when the sample 
is small in comparison with the length of the electromag- 
netic wave. 

The imaginary part of the dielectric susceptibility e" 
affects the Q-factor and results in a diminishing of the 
microwave power U transmitted through the cavity. The 
relative change of U is given by the relation: 



(3) 



The analogous relations could be used to determine the 
magnetic permeability, ji = fjf + if/', in the case when 
the sample is placed at a maximum of the microwave 
magnetic field. 

The sample used for the microwave susceptibility mea- 
surements had the dimensions 1x1x0.3 mm and, thus, 
was much smaller than the half of the length of the stand- 
ing electromagnetic wave in the cavity (5mm). This fact 
enables us to separate the magnetic and dielectric sus- 
ceptibilities by positioning the sample at different points 
of the cavity. The a- and 6-directions were lying in the 
plane of the plate, the c-axis being perpendicular to the 
plate. 

The change of the refractive index n with temperature 
was measured in the frequency range between 1.2 and 3.0 
THz by registrating the interference pattern at different 
temperatures in the transmittance spectra of thin (14- 
110 mkm) plates of NaV20s single crystals cleaved per- 
pendicular to the c-direction. The spectra were measured 
using the BOMEM DA3.002 Fourier transform spectrom- 
eter at the resolution of 0.2-2.0 cm -1 . The incident light 
was polarized either along the a— axis (E || a) or along 
the 6-axis (E || b). We have also performed E || c mea- 
surements at room temperature. 

The position v m (in wavenumbers) of the m-th maxi- 
mum in the transmittance spectrum is given by the ex- 
pression 



m 

2d^' 



(4) 



where d is the thickness of the plate. We determined the 
three principal values of the refractive index along the 
a-, b- 1 and c-axes by measuring the distance between the 
adjacent maxima of the interference pattern. 

The relative changes of the dielectric constant e' were 
found from the measured shifts of the interference max- 
ima according to the relation 



— - 2— - -2— 

e' n v m 



(•5) 



The dielectric constant is the total polarizability of the 
sample divided by the sample volume. The change of 
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the sample size at the phase transition was observed by 
studying the thermal expansion.E3 The relative change of 
the sample size in the vicinity of the phase transition was 
of the order of 10~ 4 . In our microwave measurements the 
change of the resonator frequency is caused by the total 
polarizability of the sample. Therefore the frequency of 
the resonator is not affected by the change of the sample 
size if it takes place without the change of the total polar- 
izability. Because we calculate the value of Se' supposing 
the constant size of the sample, the experimental data 
on the microwave susceptibility can not be attributed to 
the thermal expansion. 

On the contrary, the change of the refractive index 
measured at far infrared frequency by the observation of 
the interference pattern should be influenced by a change 
of the sample thickness and volume. In the case of the 
constant total polarizability the resulting relative change 
of the refractive index is to be of the order of the relative 
change of the size. 

Single crystals of stoichiometric NaX^Os used in our 
experiments have been obtained by a melt growth 
method using NaVC>3 as a flux. The details of the growth 
procedure are described in Ref. [2^. Our samples were of 
the same growth procedure as used in Ref. 2fl The tem- 
perature T c of the phase transition obtained from the be- 
ginning of the magnetic susceptibility drop is 36±0.5K. 



III. EXPERIMENTAL RESULTS 
A. Microwave measurements 

The shift of the resonant frequency of the cavity pro- 
duced by the sample at T=40K was of about 400 MHz 
both for ~E mw || b and E mji , || a. It corresponds to the 
dielectric constant values e' ab =12 ±1. 

For E mu , || c, when the microwave electric field is per- 
pendicular to the plate, the shift of the resonator fre- 
quency was only 40 ±15 MHz. Either the small value 
of e' c or the depolarization effect could result in this fre- 
quency shift which is much smaller than in the case of 
E mu , || a or E mw \\ b. Thus, the absolute value of e' c 
could not be estimated accurately enough: following the 
formula (2) we obtain the value e' c between 4 and infin- 
ity. However, the frequency shift of the resonator and the 
value of relative change of e c at the phase transition can 
be measured quite accurately. We used the value e' c =8 
obtained in quasistatic measurementstU and in the far in- 
frared experiments described below, for the calculations 
of Se c according to the formula (2). 

The temperature dependence of the change of the di- 
electric function was derived from the dependence of the 
resonator frequency shift with respect to the frequency 
at T=1.5 K. The value of the change of the imaginary 
part was derived from the temperature dependence of 
the transmitted signal relative to the reference value of 
this signal at 40 K. 



The temperature dependence of the change of the di- 
electric susceptibility is shown in Fig. 1. The real part 
of the dielectric constants e' a c decreases below the phase 
transition temperature, while the real part of eb does not 
show any observable change. The imaginary part of ex- 
hibits a well-pronounced peak. The change of the imagi- 
nary part of e a and e c were found to be smaller than the 
apparatus noise level. 

A change in the response of the sample to the mi- 
crowave magnetic field was only found for the imagi- 
nary part of the magnetic susceptibility. These results 
are shown in Fig. 2. 

We reexamined the CuGe03 crystals studied in 
Refs. |2^,|23| and have not found any changes of the dielec- 
tric and magnetic susceptibilities of comparable values. 



B. Far infrared measurements 

For the frequency of 1.2 THz we have obtained the 
following values of e' at room temperature: e^=15.0±0.6, 
e' b =10.2±0.2, e' c =7.5±0.2. These data are close to those 
for e' at zero frequency found earlier from the fitting of 
the reflectivity spectra of NaV2CU-,single crystals by the 
model of independent oscillators.E-3 

Fig. 3 shows the temperature dependence of the rela- 
tive change of the refractive index at different far infrared 
frequencies for the electric field polarized along the a- and 
&-axes. Simultaneously with diminishing of the refractive 
index n a , the absorption present in the low-frequency re- 
gion of the far infrared spectrum of NaV20s above T c 
decreases markedly when temperature decreases below 
T c (see the inset of Fig. 3 and Ref. [l6|). We checked that 
despite these changes of the absorption coefficient, equa- 
tion (5) allows us to determine the changes of n a with 
sufficiently good precision. The relative change of the 
refractive index is of about 0.05 and corresponds to the 
change of the dielectric constant that is approximately 
twice as large as the value found at the microwave fre- 
quency of 36 GHz. 

We have not found any observable changes of e' b (or 
rib). Unfortunately, we could not observe the interference 
pattern for E || c polarization at low temperatures and, 
hence, to detect the possible changes of e' c . 



IV. DISCUSSION 

The observed dielectric anomalies could, in princi- 
ple, arise from—G) an anomalous thermal expansion of 
NaV205 at T C J13 (ii) a change of the lattice polarizabil- 
ity and (iii) a rearrangement of electrons below T c . 

As it was described in Sec. II the microwave measure- 
ments of e are insensitive to changes of a sample size, thus 
the reason (i) could not result in the dielectric constant 
data shown in_Fig. 1. Besides that, the thermal expan- 
sion anomalycJ could result in changes of e' of the order 
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of 10~ 4 which is 2 orders of magnitude smaller than the 
observed value both at microwave and far infrared fre- 
quencies. 

As for (ii), the changes in the lattice polarizability 
should be accompanied by the corresponding changes in 
the phonon spectrum below T C S3 Our estimates show 
that the new modes that appear below T c cannot ac- 
count for the observed change of e' due to their relatively 
small oscillator strengths. 

The observed dielectric anomalies of relatively large 
amplitudes agree with the conclusion about the charge 
ordering at T C =35K in NaV2C>5. Analyzing the form of 
the temperature dependence of e one can choose an ad- 
equate model of the charge ordering. The model, with a 
localization of electrons on one leg of the ladderEI corre- 
sponds to the ferroelectric ordering with the spontaneous 
electric polarization along the a— axis. In this case the 
peak in the temperature dependence of e' a should be ob- 
served. Since the anomalies found in the real part of 
the dielectric constant have the form of steps (typical 
for the antiferroelectric transition) rather than peaks, we 
conclude that the charge ordering occurs without a gen- 
eration of a macroscopic dipole moment. The charge or- 
dering of an antiferroelectric type may be in the form 
of the zigzag distribution of V 4+ -ions as described in 
Refs. In this case the decrease of the dielectric 

constant e' a can be attributed to the shrinking of the 
electron clouds along the a— direction due to the local- 
ization of electrons on the distinct V-sites, instead of 
being smeared along the V-O-V orbitals. The observed 
diminishing of the optical absorption coefficient in the 
frequency range studied (1.2-3.0 THz) below T c is prob- 
ably also connected with the above mentioned shrinking 
of the electron clouds, the absorption itself being of an 
electronic origin (see, in particular, Rcf. p4| ). 

The coupled spin and charge degrees of freedom in 
vanadium ladders were described in Ref. ^ using the 
spin-isospin Hamiltonian. The isospin concept is intro- 
duced for the description of a charge ordering: Isospin 
t[- = +1/2 corresponds to the electron localized on one 
end of the rung i, while r? = —1/2 describes the elec- 
tron localized on the other end. A charge ordering would 
result in the formation of the dipole moments d\ along 
a-direction on the rungs, with d\ oc rf . The zigzag or- 
dering corresponds here to r z -antiferromagnetism. The 
coupling of an external electric field along the a-axis to 
the local (rung) order parameter is given by the term 
—E z t*, similar to the coupling —H z sf in magnetic sys- 
tems, and the resulting behavior of e' a should be similar 
to the behavior of the parallel magnetic susceptibility, x\\ j 
in an antiferromagnet. This consideration is in a qual- 
itative agreement with the experimental data shown in 
Fig. 1. 

Using the Landau mean field approach near T c , one 
can show that the decrease of the dielectric susceptibility 
is proportional to the square of the charge order param- 
eter, i.e., the value of Se' a is proportional to the intensity 
of the additional x-ray reflexes that appear at T c . As 



shown in the inset of Fig. 1, there is, indeed, a good cor- 
relation below ~ (T c — 2X)between 5e' a and the x-ray 
intensity measured in Ref. Eq. There is some deviation 
of the data close to and above T c - the value of 5e' a does 
not vanish at T c and shows a tail up to T — T c w 8K. 
This deviation above T c can naturally be ascribed to the 
existence of short-range charge ordered areas which are 
fluctuating in some frequency range. These fluctuations 
give rise to the dynamic dielectric susceptibility above T c , 
but are averaged out and vanish when the susceptibility 
is measured on a large time scale as it was observed at 
lKHz in Ref. [D§ 

The nature of the behavior of the dielectric constant 
e' c is less clear, especially, in view of the difference in the 
anomaly of e' c obtained in our measurements (see Fig. 1) 
and in the quasistatic measurements of Ref. where 
the A— type anomaly was observed for e' c instead of the 
step- like behavior shown in Fig. 1. 

Both these measurements manifest the appearance of 
local dipole moments along the c-direction at the transi- 
tion temperature. Possible origin of such moments may 
be, e.g., the shifts of the V-ions within 05-pyramides 
along c-direction. The charge ordering (localization of 
the electrons on particular V ions) accompanied by the 
shifts of these ions along the c-direction can modify the 
local dipole moments in c-direction. 

From the crystal structure of NaV205 one would then 
again expect an antiferroelectric ordering of such dipoles, 
which would give the behavior of e' c consistent with the 
one observed in the present work (see Fig. 1). The reason 
for different type of the e' c anomaly observed in our mea- 
surements and in quasistatic measurements of Ref. [l8| 
is not clear at present. A possible explanation might 
be related to the small but measurable conductivity of 
NaV2C>5 single crystals wjhich has a peak of nearly 50 per- 
cents in amplitude at T c .c3 The conductivity contributes 
only to the imaginary part of the microwave dielectric 
constant and at 40 K the contribution is of the order 
of 10~ 7 . Thus, the conductivity is too small to affect 
the results of the microwave measurements. However, 
the influence of this small conductivity on the, results of 
the electric capacity measurements at lKHzlia could be 
significant, because the observed value of conductivity 
results in the capacitor discharge time of about 10 -3 sec. 
Therefore, the electric field potential within the capac- 
itor might be disturbed by a leakage current and this 
fact might result in an additional change of the capacity. 
Further experiments are needed to elucidate this point 
and also to check whether the anomalies in e' c are indeed 
related to the shifts of V ions along the c-axis below T c 
as argued above. 

The picture of the phase transition in NaV20s as be- 
ing mainly of the charge-ordering type with some degree 
of a short-range order persisting above T c is also qualita- 
tively consistent with the behavior of the dielectric losses 
e" shown in Fig. 1. There should definitely exist slow 
enough charge fluctuations close to T c , which would give 
rise to the absorption of the microwave power. Why is 
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the anomaly in e" the most pronounced for the electric 
field directed along b- axis is not clear yet. 

Correspondingly, a coupling of the charge (and lat- 
tice) degrees of freedom to spins, definitely present in 
this compound, should result in magnetic losses, which 
were indeed observed (see Fig. 2). This coupling orig- 
inates from the Pauli principle, which relates the sym- 
metry of the spatial and spin parts of the wavefunction 
of two exchange-coupled electrons from the neighboring 
rungs of the ladder .tij As a result, triplet excitations on 
the ladder rungs could carry a dipole moment x3 

Note that the peak in the imaginary part of the mag- 
netic susceptibility is of significant value and exceeds sev- 
eral times the value of the real part of the susceptibility 
at T c . Therefore, the fluctuations of the spin gap alone 
could not explain the peak of magnetic losses and such 
peak is indeed not observed in CuGe03, where the charge 
ordering is absent. 



V. CONCLUSIONS 

The anomalies of the dielectric and magnetic suscepti- 
bility in the microwave and far infrared frequency range 
were found. These anomalies are in agreement with the 
zigzag model of the charge ordering at the phase tran- 
sition in NaV2C>5. The observation of the pronounced 
dielectric constant anomaly and of the peak in magnetic 
losses at T c indicates that the nature of the spin gap 
opening in NaV 2 5 is different from that of CuGe0 3 . 
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Figure captions 

Fig.l Temperature dependence of the real and imag- 
inary parts of the dielectric constant at the frequency 
36 GHz. The inset demonstrates the comparison of the 
measurements of Se' a (present work) and of the intensity 
/ of the x-ray reflexes^: a = 1 - 5e' a (T)/5e' a (45K). 

Fig. 2 Temperature dependence of the change of the 
imaginary part of the magnetic permeability at the fre- 
quency 36 GHz. 

Fig. 3 Temperature dependence of the relative change 
of the refractive index n a at the frequency 1.1 THz (□), 
1.4 THz (A), 1.9 THz (V), 2.2 THz (O), 2.6 THz (o), 
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2.8 THz (x), 3.1 THz (+) and n b at 3.1 THz (•). Inset 
presents the temperature dependence of the integrated 

12 THz 

intensity of the far infrared absorption I = J (3{y)dv, 

1.0 THz 

(3 being the absorption coefficient. 
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Fig. 1 Smirnov et.al 
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Fig. 2 Smirnov et al 




Fig. 3 Smirnov et. al. 



